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Abstract 

Let B he a. centrally symmetric convex polygon of and ||p — q|| be the distance 
between two points p, q € in the normed plane whose unit ball is B. For a set T of n 
points (terminals) in R^, a B-Manhattan network on T is a network N{T) = {V, E) with 
the property that its edges are parallel to the directions of B and for every pair of terminals 
ti and tj, the network N{T) contains a shortest B-path between them, i.e., a path of 
length ||ti — tj ||. A minimum B-Manhattan network on T is a S-Manhattan network of 
minimum possible length. The problem of finding minimum i3-Manhattan networks has 
been introduced by Gudmundsson, Levcopoulos, and Narasimhan (APPROX'99) in the 
case when the unit ball B is a square (and hence the distance ||p — q|| is the li or the 
Zoo-distance between p and q) and it has been shown recently by Chin, Guo, and Sun 
[6] to be strongly NP-complete. Several approximation algorithms (with factors 8,4,3, 
and 2) for the minimum Manhattan problem are known. In this paper, we propose 
a factor 2.5 approximation algorithm for the minimum B-Manhattan network problem. 
The algorithm employs a simplified version of the strip-staircase decomposition proposed 
in our paper [5] and subsequently used in other factor 2 approximation algorithms for the 
minimum Manhattan problem. 

Keywords. Normed plane, distance, geometric network design, Manhattan network, approxi- 
mation algorithms. 



I Introduction 

1.1 Normed planes 

Given a compact, centrally symmetric, convex set B in the plane M^, one can define a norm 

I I • 1 1 := 1 1 • I |e : — 7> by setting 1 1 v| | = A, where v = Au and u is a unit vector belonging to 
the boundary of B. We can then define a metric d := d^ on by setting d[p,q) = ||p — q||- 
The resulting metric space (M^,^^) is called a normed (or Minkowski) plane with unit disk 
(gauge) B [2, 24]. In this paper, we consider normed planes in which the unit ball ;B is a 
centrally symmetric convex polygon (i.e., a zonotope) of M^. We denote by bo, . . . b2m-i the 
2m vertices of B (in counterclockwise order around the circle) as well as the 2m unit vectors 
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that define these vertices. By central symmetry of B, = — bfc_|_m for k = 0, . . . ,m — 1. A 
legal k-segment of (M^, dg) is a segment pq lying on a line parallel to the line passing via 
and hk+m- A legal path 7r(p, q) between two points p, q of M? is any path connecting p and q 
in which all edges are legal segments. The length of 7r(p, q) is the sum of lengths of its edges. 
A shortest B-path between p and q is a legal (p, q)-path of minimum length. The best known 
example of normed planes with polygonal unit balls is the /i-plane (also called the rectilinear 
plane) with norm ||v|| = |a;(v)| + |y(v)|. The unit ball of the /i-plane is a square whose two 
diagonals lie on the x-axis and y-axis, respectively. The li-distance between two points p and 
q is d{p,q) := ||p — q||i = \x{p)—x{(i)\ + \y{p) — y{q)\. The legal paths of the rectilinear plane 
are the paths consisting of horizontal and vertical segments, i.e., rectilinear paths. Another 
important particular case of polygonal norms is that of A-norms (alias uniform norms) [4, 3] 
for which the unit ball ;S is a regular polygon. 

1.2 Minimum Manhattan and i3-Manhattan network problems 

A rectilinear network N = {V, E) in consists of a finite set V of points and horizontal and 
vertical segments connecting pairs of points of V. The length of N is the sum of lengths of its 
edges. Given a finite set T of points in the plane, a Manhattan network [15] on T is a rectilinear 
network N(T) = (V, E) such that TCI/ and for every pair of points in T, the network N{T) 
contains a shortest rectilinear path between them. A minimum Manhattan network on T is 
a Manhattan network of minimum possible length and the Minimum Manhattan Network 
problem (MMN problem) is to find such a network. 

More generally, given a zonotope B cM?, a network N = (V, E) consists of a finite set 
V of points and legal segments connecting pairs of points of V (the edges of N). The length 
1{N) of is the sum of lengths of its edges. Given a set T = {ti, . . . , t„} of ?i points (called 
terminals), a B-Manhattan network on T is a ;S-network A^(r) = iy, E) such that T 'ZV and 
for every pair of terminals in T, the network N{T) contains a shortest ;S-path between them 
(see Fig. 1). A minimum B-Manhattan network on T is a S-Manhattan network of minimum 
possible length and the Minimum ;S-Manhattan Network problem {B- MMN problem) is to find 
such a network. Fig. 3 illustrates the evolution of a minimum ;B-Manhattan network defined 
on the same set of terminals when the number of directions in the unit ball B is increasing 
(the directions of the unit ball are indicated at the upper left corner of each figure). 

1.3 Known results 

The minimum Manhattan network problem has been introduced by Gudmundsson, Levcopou- 
los, and Narasimhan [15]. Gudmundsson et al. [15] proposed an 0(n^)-time 4-approximation 
algorithm, and an 0(n log ?i)-time 8-approximation algorithm. They also conjectured that 
there exists a 2-approximation algorithm for this problem and asked if this problem is NP- 
complete. Quite recently, Chin, Guo, and Sun [6] solved this last open question from [15] and 
established that indeed the minimum Manhattan network problem is strongly NP-complete. 
Kato, Imai, and Asano [17] presented a 2-approximation algorithm, however, their correct- 
ness proof is incomplete (see [1]). Following [17], Benkert, Wolff, Shirabe, and Widmann 



2 




tion for this instance does not belong 
Figure 1: A S-Manhattan network in the normed to the grid T (the unit baU ^ is a 
plane whose unit ball is depicted in Fig. 4 hexagon) 

[1] described an 0(n log n)-time 3-approximation algorithm and presented a mixed-integer 
programming formulation of the MMN problem. Nouioua [20] and later Fuchs and Schulze 
[12] presented two simple 0(n log ?i)-time 3-approximation algorithms. The first correct 2- 
approximation algorithm (thus solving the first open question from [15]) was presented by 
Chepoi, Nouioua, and Vaxes [5]. The algorithm is based on a strip-staircase decomposition 
of the problem and uses a rounding method applied to the optimal solution of the flow based 
linear program described in [1]. In his PhD thesis, Nouioua [20] described a 0(n log n)-time 
2-approximation algorithm based on the primal-dual method from linear programming and 
the strip-staircase decomposition. In 2008, Guo, Sun, and Zhu [13, 14] presented two com- 
binatorial factor 2 approximation algorithms, one with complexity 0{'n?) and another with 
complexity O(rilogn) (see also the PhD thesis [22] of Schulze for yet another 0(n log n)-time 
2-approximation algorithm). Finally, Seibert and Unger [21] announced a 1.5-approximation 
algorithm, however the conference format of their paper does not permit to understand the 
description of the algorithm and to check its claimed performance guarantee (a counterexam- 
ple that an important intermediate step of their algorithm is incorrect was given in [12, 22]). 
Quite surprisingly, despite a considerable number of prior work on minimum Manhattan net- 
work problem, no previous paper, to our knowledge, consider its generalization to normed 
planes. 

Gudmundsson et al. [15] introduced the minimum Manhattan networks in connection 
with the construction of sparse geometric spanners. Given a set T of n points in a normed 
plane and a real number t > 1, a geometric network is a t-spanner for T if for each pair of 
points p, q G T, there exists a (p, q)-path in N of length at most t times the distance ||p — q| 
between p and q. In the Euclidian plane and more generally, for normed planes with round 
balls, the line segment is the unique shortest path between two endpoints, and therefore the 
unique 1-spanner of T is the complete graph on T. On the other hand, if the unit ball of 
the norm is a polygon, the points are connected by several shortest ^B-paths, therefore the 
problem of finding the sparsest 1-spanner becomes non trivial. In this connection, minimum 
;B-Manhattan networks are precisely the optimal 1-spanners. Sparse geometric spanners have 
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Figure 3: "Morphing" a minimum ;B-Manhattan network 



applications in VLSI circuit design, network design, distributed algorithms and other areas, 
see for example the survey of [10] and the book [18]. Lam, Alexander sson, and Pachter 
[19] suggested to use minimum Manhattan networks to design efficient search spaces for pair 
hidden Markov model (PHMM) alignment algorithms. 

Algorithms for solving different distance problems in normed spaces with polygonal and 
polyhedral balls were proposed by Widmayer, Wu, and Wang [25] (for more references and 
a systematic study of such problems, see the book by Fink and Wood [11]). There is also an 
extensive bibliography on facility location problems in normed spaces with polyhedral balls, 
see for example [9, 23]. Finally, the minimum Steiner tree problem in the normed planes was 
a subject of intensive investigations, both from structural and algorithmic points of view; 
[4, 3, 8] is just a short sample of papers on the subject. 

2 Preliminaries 

2.1 Definitions, notations, auxifiary results 

We continue by setting some basic definitions, notations, and known results. Let ;B be a 
zonotope of with 2m vertices bg, . . . b2m-i having its center of symmetry at the origin 
of coordinates (see Fig. 4 for an example). The segment Sk ■= i>k^k+i(mod 2m) is a side 
of B. We will call the line ii passing via the points and hk+m an extremal line of B. 
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Figure 4: A unit ball B and an interval I(p, q) 

Two consecutive extremal lines ik and ik+i defines two opposite elementary k- cones Ck and 
Cfc+m, = —Ck containing the sides Sk and Sk+m-, respectively. We extend this terminology, 
and call elementary k-cones with apex v the cones Cfc(v) = + v and — Cfc(v) = Ck+m + v 
obtained by translating the cones Ck and Ck+m by the vector v. We will call a pair of 
consecutive lines Dk = {(-k, ^k+i} a direction of the normed plane. Denote by B{^, r) = r -S+v 
the ball of radius r centered at the point v. 

Let /(p, q) = {z G : d(p, q) = (i(p, z) + (i(z, q)} be the interval between p and q. The 
inclusion pq C /(p, q) holds for all normed spaces. If B is round, then pq = /(p,q), i.e., 
the shortest path between p and q is unique. Otherwise, /(p, q) may host a continuous set 
of shortest paths. The intervals /(p, q) in a normed plane (and, more generally, in a normed 
space) can be constructed in the following pretty way, described, for example, in the book 
[2]. If pq is a legal segment, then pq is the unique shortest path between p and q, whence 
-^(Pil) = PQ- Otherwise, set r = d(p, q). Let s'^. be the side of the ball S(p, r) containing 
the point q and let s'- be the side of the ball S(q, r) containing the point p. Notice that 
these sides are well-defined, otherwise q is a vertex of ^(p, r) and pq is a legal segment. The 
segments s'^ and s'- are parallel, thus |A; — j| = m, say k < m and j = k + m. Then /(p, q) 
is the intersection of the elementary cones Ck{p) and Ck+mi^) = —Ck{(l) (see Fig. 4 for an 
illustration): 

Lemma 2.1 [2] /(p,q) = C7fc(p) n (-Cfc(q)). 



An immediate consequence of this result is the following characterization of shortest B- 
paths between two points p and q. 

Lemma 2.2 If pq is a legal segment, then pq is the unique shortest B-path. Otherwise, 
i//(p, q) = Ck{p) n (— Cfc(q)), then any shortest B-path '/r(p, q) between p and q has only 
k-segments and {k + \)-segments as edges. Moreover, '7r(p, q) is a shortest B-path if and only 
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if it is monotone with respect to Ik and It+i^ ^.e., the intersection o/7r(p,q) with any line i 
parallel to the lines ik,ik+i is empty, a point, or a (legal) segment. 

Proof. The first statement immediately follows from Lemma 2.1. Suppose that pq is not a 
legal segment and I(p, q) = Cfc(p) fl Cfc+m,(q)- Let uv be the first edge on a shortest path 
7r(p, q) from p to q which is neither a A:-segment nor a (A; + l)-segment. Since u G /(p, q) = 
Cfe(p) n C'/c+m(q)) the point q belongs to the cone Cfc(u) and the point u belongs to the 
cone Cfc+m(q), whence I(u, q) = Cfe(u) fl Cfc-|_m(q)- Obviously, the point v belongs to I(u, q). 
However, by the choice of the segment uv and the fact that c^nd ik+i are consecutive 
lines that forms a direction, the point v cannot belong Cfc(u), a contradiction. This shows 
that any shortest legal path 7r(p, q) between p and q has only k- and {k + l)-segments as 
edges. Additionally, the intersection of 7r(p, q) with any line i parallel to ik or i^+i is empty, 
a point, or a (legal) segment. Indeed, pick any two points in this intersection. Since the 
legal segment defined by these points is the unique shortest path between them, it must also 
belong to the intersection of 7r(p, q) with i. Conversely, consider a monotone path 7r(p, q) 
between p and q, namely suppose that the intersection of 7r(p, q) with any line i parallel to 
the lines ik or i^-^-i is empty, a point, or a (legal) segment). We proceed by induction on 
the number of edges of ■7r(p, q). The monotonicity of 7r(p, q) implies that 7r(p, q) lies entirely 
in the interval /(p,q). In particular, the neighbor u of p in vr(p, q) belongs to /(p,q). The 
subpath 7r(u, q) of 7r(p, q) between u and q is monotone, therefore by induction assumption, 
7r(u, q) is a shortest path between u and q. Since pu is a legal segment and u G /(p, q), we 
immediately conclude that vr(p, q) is also a shortest path between p and q. □ 

We continue with some notions and notations about the S-MMN problem. Denote by 
OPT(r) the length of a minimum ;B-Manhattan network for a set of terminals T. For a 
direction Dk = {ik,ik+i}, denote by Fk the set of all pairs {i,j} (or pairs of terminals 
{tj,tj}) such that any shortest ;B-path between tj and tj uses only fc-segments and (fc + 1)- 
segments. Equivalently, by Lemma 2.2, consists of all pairs of terminals which belong to 
two opposite elementary cones Cfc(v) and — Cfc(v) with common apex. For each direction 
Dk and the set of pairs F^, we formulate an auxiliary problem which we call Minimum 1- 
Directional Manhattan Network problem (or l-DMMN(Fk) problem): find a network N'^^(T) 
of minimum possible length such that every edge of N'^^'{T) is an fe-segment or an [k + 1)- 
segment and any pair {tj, tj} of F^ is connected in N^^{T) by a shortest B-path. We denote 
its length by OPTfc(T). We continue by adapting to 1-DMMN the notion of a generating 
set introduced in [17] for MMN problem: a generating set for F^ is a subset F of F^ with 
the property that a ;B-Manhattan network containing shortest S-paths for all pairs in F is a 
1-Directional Manhattan network for F^.- 

2.2 Our approach 

Let N*{T) be a minimum S-Manhattan network, i.e., a ;B-Manhattan network of total length 
1{N*{T)) = OPT(r). For each direction Dk, let N*{T) be the set of fc-segments and {k + 1)- 
segments of N*{T). The network N^iT) is an admissible solution for 1-DMMN (Ffc), thus 
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the length 1{N^{T)) of N^{T) is at least OPTfc(r). Any fc-segment of N*{T) belongs to 
two one-directional networks N^{T) and N^_-^^(T). Vice- versa, if Nk{T),k = 0, . . . ,m — 1, 
are admissible solutions for the l-DMMN(i<fc) problems, since IJfc^^ -^fc = T x T, the 
network N{T) = \^]^=q Nk{T) is a ;B-Manhattan network. Moreover, if each Nk{T) is 
an a-approximation for respective 1-DMMN problem, then the network N{T) is a 2a- 
approximation for the minimum ,B-Manhattan network problem. Therefore, to obtain a 
factor 2.5-approximation for S-MMN, we need to provide a 1.25-approximation for the 1- 
DMMN problem. The remaining part of our paper describe such a combinatorial algorithm. 
The 1-DMMN problem is easier and less restricted than the ;S-MMN problem because we 
have to connect with shortest paths only the pairs of terminals of the set Ff^ corresponding 
to one direction D^, while in case of the MMN problem the set T x T of all pairs is parti- 
tioned into two sets corresponding to the two directions of the ii-plane. For our purposes, 
we will adapt the strip-staircase decomposition of [5], by considering only the strips and the 
staircases which "are oriented in direction D^". 

3 One-directional strips and staircases 

In the next two sections, we assume that = {Z^, Z^+i} is a fixed but arbitrary direction of 
the normed plane. We recall the definitions of vertical and horizontal strips and staircases 
introduced in [5] . Then we consider only those of them which correspond to pairs of terminals 
from the set F^, which we call one-directional strips and staircases. We formulate several 
properties of one-directional strips and staircases and we prove those of them which do not 
hold for usual strips and staircases. 

Denote by and Lfc+i the set of all lines passing via the terminals of T and parallel to 
the extremal lines £k and £k+i, respectively. Let be the grid defined by the lines of Lk and 
Lfc+i. The following lemma can be proved in the same way as for rectilinear Steiner trees or 
Manhattan networks (quite surprisingly, this is not longer true for the ;6-MMN problem: Fig. 
2 presents an instance of jB-MMN for which the unique optimal solution does not belong to 
the gridF:=Ur=o'rA:): 

Lemma 3.1 There exists a minimum 1-Directional Manhattan Network for contained in 
the grid F^ . 

For two terminals ti,tj/, set Ri^ii := /(tj,tj/). A pair tj,tj/ defines a k-strip Ri i' if either 
(i) (degenerated strip) tj and tj/ are consecutive terminals belonging to the same line of 
or (ii) ti and tj/ belong to two consecutive lines of and the intersection of Rn' with any 
degenerated fc-strip is either empty or one of the terminals tj or tj/; see Fig. 6 of [5]The two 
/c-segments of i?j^j/ are called the sides of -Rj^j/. The (k + 1) -strips and their sides are defined 
analogously (with respect to Lfc+i)- With some abuse of language, we will call the fc-strips 
horizontal and the {k + l)-strips vertical. If a pair {tj,tj/} defining a horizontal or a vertical 
strip Ri^i' belongs to the set Fk, then we say that i?j.j/ is a one- directional strip or a 1-strip, 
for short. Denote by Fj^ the set of all pairs of F^ defining one-directional strips. 
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Figure 5: Strips, staircases, and completion 



Lemma 3.2 If Ri^i' and Rjji are two horizontal 1-strips or two vertical 1-strips, then Ri^ii PI 
^ij' = ^if n {j,/} = and Ri,y n i?jjv = {ti} if {i, i'] n {j,/} = {i}. 

Proof. From the definition follows that if Ri^ii and Rjji are both degenerated or one is 
degenerated and another one not, then they are either disjoint or intersect in a single terminal. 
If Ri^i' and -Rjj' are both non-degenerated and intersect, then from the definition immediately 
follows that the intersection is one point or a segment belonging to their sides. However, if 
Ri^ii and Rjji intersects in a segment, then one can easily see that at least one of Ri^ii and 
Rj^ji cannot be a 1-strip. □ 

We say that a vertical 1-strip Ri ii and a horizontal 1-strip Rj ji (degenerated or not) form 
a crossing configuration if they intersect (and therefore cross each other). 



Lemma 3.3 If Ri^ii and Rjj' form a crossing configuration, then from the shortest B-paths 
between ti and tj/ and between tj and tj/ one can derive shortest B-paths connecting tj,tj/ 
and tif ,tj, respectively. 

For a crossing configuration defined by the 1-strips Ri^i' , Rjj' , denote by o and o' the two 
opposite corners of the parallelogram i?j j/ PI Rjj', such that the cones Cfc(o) and — Cfc(o') do 
not intersect the interiors of Ri^i' and Rjj'- Additionally, suppose without loss of generality. 
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that tj and tj belong to the cone Ck{o), while tj^ and tji belong to the cone — Cfc(o'). Denote 
by Tij the set of all terminals G (T\{ti, tj})nCfc(o) such that (— Cfc(t;))\(— Cfc(o)) does not 
contain any terminal except ti. Denote by Sij^i/j/ the region of Cfc(o) which is the union of the 
intervals I{ti,o),ti € Tjj, and call this polygon an one- directional staircase or a 1-staircase, 
for short; see Fig. 5 and Figures 7,8 of [5] for an illustration. Note that Sij^i/ji is bounded 
by the 1-strips Ri^ii and Rjj' and a legal path between tj and tj passing via all terminals 
of Tij and consisting of A;-segments and {k + l)-segments. The point o is called the origin 
and Ri^ii and Rjj' are called the basis of this staircase. Since I(t;, o) C {—Ck{ti)) \ (— Cfc(o)) 
for all ti € Tj j, /(t;,o) fl T = {t;} and therefore 5jj|j/jv n T = Tij. For the same reason, 
there are no terminals of T located in the regions Q' and Q" depicted in Fig. 5 (Q' is the 
region comprised between the leftmost side of i?j^ j/ , the highest side of Rjj' , and the line of 
Lfc passing via the highest terminal of Tjj, while Q" is the region comprised between the 
rightmost side of i?j,j', the lowest side of Rjj', and the line of Lfc+i passing via the rightmost 
terminal of Tjj). Analogously one can define the set Tj/j/ and the staircase Si/j/^ij with 
origin o' and basis Riy and Rjj'. 

Lemma 3.4 // a 1-strip Ri^ii intersects a 1-staircase Si/j/^ij and Ri^i' is different from the 
1-strips Ri^i' and Rjj', then Ri^i' n5j/ .,/|j j is a single terminal. 

Proof. If a 1-strip Ri i' traverses a staircase Si/j/^ij, then one of the terminals t/,!^/ must 
be located in one of the regions Q' and Q", which is impossible because {Q' U Q") fl T = 0. 
Thus, if Ri^ii and 5j/ j/|j j intersect more than in one point, then they intersect in a segment 
s which belongs to one side of Ri^i' and to the boundary of Si'ji^ij. If say the 1-strip Ri^i' is 
horizontal, then necessarily s is a part of the lowest side of Ri^i' and of the highest horizontal 
side of Si'j'^ij. Let t be the highest terminal of Tij. Then either t belongs to Ri^i' and is 
different from ti,ti', contrary to the assumption that R^i' is a strip, or t together with the 
lowest terminal t// of Ri^i' define a degenerated strip with t// belonging to Q' , contrary to the 
assumption that Q' n T = 0. □ 

Lemma 3.5 Two 1-staircases either are disjoint or intersect only in common terminals. 

Proof. From the definition of a staircase follows that the interiors of two staircases are 
disjoint (for a short formal proof of this see [5]). Therefore two staircases may intersect only 
on the boundary. In this case, the intersection is either a subset of terminals of both staircases 
or a single edge. In the second case, one of the two staircases necessarily is not a 1-staircase 
with respect to the chosen direction. □ 

Let Fj! be the set of all pairs {tj/,t;} such that there exists a 1-staircase 5jj|j/j/ with t/ 
belonging to the set Tij. The proof of the following essential result is identical to the proof 
of Lemma 3.2 of [5] and therefore is omitted. 

Lemma 3.6 F := Fj^U F^' is a generating set for F^. 
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4 The algorithm 



We continue with the description of our factor 1.25 approximation algorithm for 1-DMMN 
problem. Let Fj^ and FJ^ denote the pairs of Fj^ defining horizontal and vertical 1-strips, 
respectively. Let and ^2 be the networks consisting of lower sides and respectively upper 
sides of the horizontal 1-strips of Fj}. Analogously, let and S2 be the networks consisting 
of rightmost sides and respectively leftmost sides of the vertical 1-strips of F^ . The algorithm 
completes optimally each of the networks , S2 , Sf, and S2 , and from the set of four com- 
pletions N^,N2,N^,N2, the algorithm returns the shortest one, which we will denote by 
Nk{T) (in this respect, our algorithm has some similarity with the approach of Benkert et al. 
[1]). We will describe now the optimal completion A^^ for the network Si, the three other 
networks are completed in the same way (up to symmetry). 

An optimal completion of Si is a subnetwork Ni of extending Si {Si C Ni) of smallest 
total length such that any pair of terminals of F can be connected in A^^^ by a shortest path. 
By Lemma 3.6, to solve the completion problem for Si, it suffices to (i) select a shortest path 
7r(tj,tj/) of Ffe between each pair tj,tj/ defining a vertical 1-strip Ri^i', (ii) for each horizontal 
1-strip Rjj' find a shortest path 7r(tj,tj/) between tj and tji subject to the condition that 
the lowest side s'j j, of Rjj' is already available, (iii) for each staircase Sij^i/j/ whose sides 
are Ri^i' and Rjj' select shortest paths from the terminals of Tjj to the terminal tj' subject 
to the condition that the lowest side s'- ■, of Rj is already available. We need to minimize 
the total length of the resulting network Ni over all vertical 1-strips and all 1-staircases. To 
solve the issue (ii) for a horizontal 1-strip Rjj', we consider the rightmost 1-staircase Sij^i/ji 
having Rj j' as a basis, set Tij := Tij U {tj}, and solve for this staircase the issue (iii) for 
the extended set of terminals. For all other 1-staircases Sij^iiji and Sj_i ji\i j having Rj ji as a 
basis, we will solve only the issue (iii) for Tjj and Tii^i, respectively. 

To deal with (iii) , for each vertical 1-strip j/ , we pick each shortest path vr of Tk between 
tj and tj/, include it in the current completion, and solve (iii) for all 1-staircases having Ri^ii 
as a vertical base and taking into account that vr is already present. We have to connect the 
terminals of Tjj- by shortest paths of Vk of least total length to the terminal tj/ subject to 
the condition that the union vr U s'- is already available; see Fig. 5. For a fixed path vr, 
this task can be done by dynamic programming in 0(|Tjjp) time. For this, notice that in 
an optimal solution (a) either the highest terminal of Tjj- is connected by a vertical segment 
to s'j -I, or (b) the lowest terminal of Tjj- is connected by a horizontal segment to vr, or (c) 
Tjj- contains two consecutive (in the staircase) terminals ti,t/+i, such that t/ is connected 
to vr by a horizontal segment and t^+i is connected to s'- by a vertical segment. In each 
of the three cases and subsequent recursive calls, we are lead to solve subproblems of the 
following type: given a set T' of consecutive terminals of Tij, the path vr and a horizontal 
segment s' , connect to tji the terminals of T' by shortest paths of least total length if the 
union vr U s' is available. We define by C^^^, the optimal completion obtained by solving by 
dynamic programming those problems for all staircases having j/ as a vertical basis (note 
that vr C CJ^, however Sj^dC^^, = 0). For each vertical 1-strip Ri^i', the completion algorithm 
returns the partial completion C°^f of least total length, i.e, C°^* is the smallest completion of 
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the form C^-, taken over all 0{n) shortest paths vr running between tj and tj/ in F^. Finally, 
let be the union of all C°^* over all vertical 1-strips Ri^i' and 5"^. The pseudocode of the 
completion algorithm is presented below (the total complexity of this algorithm is 0{n^)). 

Algorithm 1 Optimal completion(5(') 

1: ^ S'^ 

2: for each vertical 1-strip Ri^i' do 

3: for each shortest path tt of F^ connecting the terminals tj and tj/ do 

4: compute the partial completion C^-, in the following way: 

6: for each 1-staircase Sij^i'j/ and each 1-staircase Si'j/^ij do 

7: if Sij^i/ji is the rightmost staircase having Rjj' as a base, then set Tj j Tij[j{tj} 

8: compute by dynamic programming the subset C of edges of F^ of least total length 

such that C U (vr U s'j j,) contains a shortest path of F^ from each terminal of Tij 
to tji or from each terminal of Tj/j/ to tj 
9: Cf,, ^ Cf,, U C 

10: end for 
11: end for 

12: let C°^f be the partial completion of least total length, i.e, C°^,* is the smallest comple- 
tion C^-, over all shortest paths vr between tj and tj/ 
13: iVf ^ U C°f* 
14: end for 
15: return NI^ 



Lemma 4.1 The network returned by the algorithm Optimal completion is an optimal 
completion for Si- 

Proof. We described above how to compute for each 1-staircase Sij^iiji a subset C of edges 
of Ffc of minimum total length such that CU (vrUs^- j,) contains a shortest path of F^ from each 
terminal of Tjj to tj/. This standard dynamical programming approach explores all possible 
solutions and therefore achieves optimality for this problem. Next, we assert that, for each 
vertical 1-strip -Rj^j/, the subset of edges C°^f computed by our algorithm, is an optimal 
completion of Si for the strip i?j^j/ and the staircases having i?j^j/ as vertical bases. Indeed, 
our algorithm considers every possible shortest path vr of F^ between tj and tj/. Once the 
path TT is fixed, the subproblems related to distinct staircases become independent and can 
be solved optimally by dynamic programming. The problems arising from distinct vertical 
1-strips are also disjoint and independent (according to Lemmas 3.4 and 3.5). Therefore the 
solution A^^* obtained by combining the optimal solutions C"^f of every vertical 1-strip i?j j/ 
is an optimal completion of S^. 

It remains to show that to obtain a completion satisfying the conditions (i),(ii), and (iii), 
it suffices for each horizontal 1-strip Rjji to add tj to the set Tj j of terminals of the rightmost 
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staircase Sij^i'j' having Rjj' as a basis and to solve (iii) for this extended set of terminals. 
Indeed, in any completion any shortest path between tj and tj' necessarily makes a vertical 
switch either before arriving at the origin o of Sij^i/ji or this path traverses the vertical basis 
of this staircase. Since the completion contains a shortest path connecting the terminals of 
the vertical basis of S'j combining these two paths, we can derive a shortest path between 
tj and t'j which turns in iZj j/ n Rjj'- As a result, we conclude that at least one shortest path 
between tj and tj' passes via o'. This shows that indeed it suffices to take into account the 
condition (ii) only for each rightmost staircase. □ 

Lemma 4.2 The network Nk{T) is an admissible solution for the problem l-DMMN{Fi?j. 

Proof. By Lemma 4.1, is a completion of Si and thus contains a shortest path between 
every pairs of vertices from F. By symmetry, we get the same result for and 
Since Nk{T) is one of these networks, by Lemma 3.6, it is admissible solution for the problem 
l-DMMN(Ffc). □ 

5 Approximation ratio and complexity 

In this section, we will prove the following main result: 

Theorem 5.1 The network Nk{T) is a factor 1.25 approximation for l-DMMN{Fk) problem 
for k = 0,...,m-l. The network N{T) := Ur=o ^fc(^) 

is a factor 2. 5 approximation for 
the B - MM N problem and can be constructed in 0{mn^) time. 

Proof of Theorem 5.1. First we prove the first assertion of the theorem. Let A/j = l{Si) = 
1{S!^) and A„ = 1{SI) = 1{S^). Further, we suppose that A;, < A^. Assume N^^ be an 
optimal 1-restricted Manhattan network for Fk- Let M be a subnetwork of -/V^^* H {Si U ^2 ) 
of minimum total length which completed with some vertical edges of N^^^ contains a shortest 
path between each pair of terminals defining a horizontal 1-strip of F^. Such M exists because 
the network N^'^ n {Si U S2) already satisfies this requirement. Further, we assume that 

l{MnS'l) > l{Mns!^). 

Lemma 5.2 1{M) = Ah- 

Proof. By Lemma 3.2, two horizontal 1-strips either are disjoint or intersect only in common 
terminals, thus any horizontal 1-strip Ri ii contributes to M separately from other horizontal 
1-strips. Since the terminals and t- defining j/ are connected in -/V^^* by a shortest path 
consisting of two horizontal segments of total length equal to the length of a side of Ri^ii and 
a vertical switch between these segments, from the optimality choice of M we conclude that 
the contribution of Ri^ii to M is precisely the length of one of its sides. □ 

Lemma 5.3 1{M D S^) > 0.5 1{M). 
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Proof. The proof fohows from the assumption 1{M fl S^) > 1{M n S2) and the fact that 
Mr\S\ and M n S"^ form a partition of M. □ 

Lemma 5.4 1{S^ \ M) < 0.25 /(iV^!^*). 

Proof. Since 1{S'^ \ M) = 1{S\) - 1{M n S^), by Lemma 5.2 and 5.3 we get l{Sl \ M) < 
0.5 1{M) = 0.5 Kh < 0.25 1{N^*). The last inequality follows from K^^*) > J^h + K (a 
consequence of Lemma 3.2) and the assumption h.^ < A„. □ 

Now, we complete the proof of Theorem 5.1. Note that 

KNk(T)) < l{Ni') = U Ni') < 1{S'^ U iV°P*) = 1{S'; \ N°f) + /(iV°P*) < 1.25 1{N°p'). 

The first inequality follows from the choice of Ni.(T) as the shortest network among the four 
completions N^,N2,N^,N2. The second inequality follows from Lemma 4.1 and the fact 
that iV°^* (and therefore 5^^ U iV°^*) is an admissible solution for the completion problem for 
Si - Finally, the last inequality follows from Lemma 5.4 by noticing that M C N^^* and thus 
KSi \ A^r*) ^ KSi \ M) < 0.25 1{N^*'). This concludes the proof of the first assertion of 
Theorem 5.1. 

Let N*(T) be a minimum ;B-Manhattan network. For each direction Dk, let N^{T) be 
the set of fc-segments and (A; + l)-segments of N*(T). By Lemma 2.2 the network N^(T) is an 
admissible solution for l-DMMN(Ffc) problem, thus 1{N^{T)) > OPTfe(r). Any A;-segment 
of N*{T) belongs to exactly two one-directional networks N^{T) and N^_^{T), we conclude 
that Er=o^OPTfc(r) < T.T=dKNk{T)) < 2 1{N*{T)) = 2 0VT{T). The first assertion of 
Theorem 5.1 implies that l{Nk{T)) < 1.25 l{Nl^^) = 1.25 OPTfc(r) for all A; = 0, . . . , m - 1, 
hence 

m—l m—l 

KN{T)) < l{Nk{T)) < 1.25 J2 OPTfc(r) < 2.5 OPT(r). 

k=0 k=0 

This concludes the proof that the approximation factor of N(T) := IJfc^^ ^k(T) is 2.5. 

To finish the proof of Theorem 5.1, it remains to analyze the complexity of the algorithm. 
First, we use a straightforward analysis to establish a O(mn^) bound on its running time. 
Then we show that this bound can be reduced to O(mn^) by using a more advanced imple- 
mentation. The time complexity of Optimal completion(S'^) is dominated by the execution of 
the dynamic programming algorithm that computes an optimal completion for each staircase 
Sij\ii^ji. The staircase Sij^i/j/ is processed 0{\Tij\) times (once for each shortest (tj, tj/)-path 
in Ffc) using a 0(|rjj|^)-time dynamic programming algorithm (each of the 0(|Tjjp) entries 
of the dynamic programming table is computed in time 0(|Tjj|)). Therefore, each stair- 
case Sij\i>j> contributes 0(|Tjj|^) to the execution of the algorithm Optimal completion(S'J*). 
Since each terminal belongs to at most two staircases, the overall complexity of Optimal 
completion(S'j*) is O(n^). This algorithm is processed to compute four optimal completions 
for each of the m directions. Therefore the total complexity of our 2.5-approximation algo- 
rithm for the ;S-MMN problem is O(mn^). 
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The following simple idea allows to reduce the contribution of each staircase Sij^^/j/ to 
0(|Tjjp) instead of 0{\Tij\'^), leading to a total complexity of O(mn^). First, note that 
among all 0(|Tjjp) subproblems, whose optima are stored in the dynamic programming 
table, only 0{\Tij\) are affected by the choice of the vr (those are the subproblems containing 
the highest and the rightmost terminal of Tjj). Therefore, instead of solving each of 0(|Tjj p) 
subproblems 0{\Tij\) times, we solve the subproblems not affected by the choice of vr only 
once. Now, consider the number of subproblems obtained by taking into account the choice of 
vr, then it is easy to verify that the total number of subproblems encountered is not OdTijl"^) 
but only 0{\Tij\'^). Since each entry of the dynamic programming table is computed in time 
0{\Tij\), we obtain a contribution of 0{\Tij\^) for each staircase Sij\j_iji, and thus a total 
complexity of 0{mn^). □ 

6 Conclusion 

In this paper, we presented a combinatorial factor 2.5 approximation algorithm for NP-hard 
minimum Manhattan network problem in normed planes with polygonal unit balls (the B- 
MMN problem). Its complexity is 0{mn^), where n is the number of terminals and 2m is 
the number of extremal points of the unit ball B. Any S-Manhattan network N{T) can be 
decomposed into m subnetworks, one for each direction of the normed plane. Each such 
subnetwork iVfc(T) ensures the existence of shortest paths between the pairs of terminals for 
which all legal paths use only k- and {k + l)-segments. We presented a factor 1.25 O(n^) 
algorithm for computing one-directional Manhattan networks, which lead to a factor 2.5 
algorithm for minimum ,8-Manhattan network problem. One of the open questions is whether 
the one- directional Manhattan network problem is NP-complete? Another open question is 
designing a factor 2 approximation algorithm for B-MMN, thus meeting the current best 
approximation factor for the classical MMN problem. Notice that polynomial time algorithm 
for 1-DMMN problem will directly lead to a factor 2 approximation for ;S-MMN. 

Notice some similarity between the 1-DMMN problem and the oriented minASS problem 
investigated in relationship with the minimum stabbing box problem [16], alias the minimum 
arborally satisfied superset problem (minASS) [7]. In the minASS problem, given a set of 
n terminals T C M^, one need to add a minimum number of points S such that for any 
pair tj,tj G T U S, either tjtj is a horizontal or a vertical segment, or the (axis-parallel) 
rectangle iijj spanned by tj,tj contains a third point of T US'. The oriented minASS problem 
is analogous to minASS problem except that the above requirement holds only for pairs 
tj,tj G T U S" such that {i,j} G Fq, i.e., tj and tj lie in the first and the third quadrants 
of the plane with the same origin. The authors of [7] presented a polynomial primal-dual 
algorithm for oriented minASS problem, however, in contrast to ;B-MMN problem, solving 
oriented minASS problems for pairs of Fq and Fi (where Fq[J F\ = (T Li S) x (T U S)) does 
not lead to an admissible solution and thus to a constant factor approximation for minASS 
(which, as we have shown before, is the case for 1-DMMN and B-MMN problems). 
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